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I. INTRODUCTION 

Since 2003, big progress on the observations of heavy- 
light mesons has been made. When checking 2014 edi¬ 
tion of Particle Data Group (PDG) [1], we notice that 
charmed meson, charmed-strange meson, bottom meson, 
and bottom-strange meson families have become more 
and more abundant, which is due to these observed can¬ 
didates of higher radial and orbital excitations of heavy- 
light mesons. In the following years, theorists and ex¬ 
perimentalists will pay more attentions to the study of 
heavy-light mesons with higher radial and orbital quan¬ 
tum numbers, especially with the running of LHCb and 
forthcoming Bellell. 

With the charmed-strange meson family as exam¬ 
ple, we introduce the research status of higher excita¬ 
tions of heavy-light meson. There are two 15 states 
(Ds(1968) and D*(2112)) and four IP states (P>:o(2317), 
Pi,i(2460), Pisi(2536), and £>*2(2573)) established in 
PDG [1]. The observed £>*1(2700) [2], £)*i(2860) [3, 4], 
and £>*3(2860) [3, 4] stimulated theorist’s interest in 
studying the properties of 25 and 1£> states [5-7], while 
the observation of £,jj(3040) [8] made us to focus on the 
2P states [9]. The research status of charmed-strange 
mesons can be found by a mini review [10] and two re¬ 
cent systematical theoretical work [6, 7]. The theoretical 
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and experimental situation of charmed meson is similar 
to that of charmed-strange meson, which can be found 
in Ref. [11]. 

Considering the above research status of heavy-light 
meson, it is suitable time to carry out the study of F-wave 
heavy-light mesons, since these IF states will be reported 
in future experiment. The calculation of mass spectrum 
of P-wave heavy-light mesons can provides valuable in¬ 
formation to experimental search for them. Before the 
present work, there were several quark model calculation 
of mass spectrum of P-wave heavy-light mesons. For ex¬ 
ample, Ebert et al. adopted the relativistic quark model 
to get the heavy-light meson spectroscopy [12], which 
includes the IP states. In Ref. [13], a relative quark 
model including the leading order correction in ljmc,b 
was applied to study heavy-light meson masses and light 
hadronic transition rates, where this study also contains 
IP states. Recently, in Refs. [6, 11], the masses of 
IP states in charmed meson and charmed-strange me¬ 
son families were obtained through the modified Godfrey- 
Isgur (GI) model, where the screening effect is considered 
in the introduced potential. For bottom and bottom- 
strange mesons, the masses of the IP states were esti¬ 
mated by the GI model in Ref. [14]. 

Although there were quark model calculations of IP 
states of heavy-light mesons, we notice that a QCD sum 
rule (QSR) study of mass spectrum of P-wave heavy- 
light mesons is still absent at present, which inspires our 
interest in performing the calculation of QSR of mass 
spectrum of P-wave heavy-light mesons. In Refs. [15, 16] 
M. A. Shifman wrote about QSR that: 

“One failure is quite obvious: the large-spin hadrons. 
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Indeed, the latter have parametrically large sizes and a 
‘sausage-like shape’ (growing with spin) and, therefore, 
it is guite clear that the basic idea of the method - ex¬ 
trapolation from short to intermediate distances - is not 
applicable. Practically, we have to stop at S=2.” 
However, it is still worth a try to applying QSR to study 
F-wslyb heavy mesons, because a) we have used the same 
method to well study D-wave heavy mesons [17] and P- 
wave heavy baryons [18]; and b) the LHCb experiments 
have just observed D-wsve heavy mesons [3, 4], and F- 
wave heavy mesons are expected in the following exper¬ 
iments. Hence, the present pioneering study not only 
provides important hint to experimental exploration of 
F-wave heavy-light mesons, but also test the applicabil¬ 
ity of QSR when applying QSR to study so higher radial 
excitations. This can be useful for quantifying potential 
overextensions of QSR in order to inspire ideas for its im¬ 
provement, especially with future experimental data on 
F-wave heavy mesons. 

The F-wave Qs {Q = c, b) heavy mesons have large 
orbital excitations L = 3, and may be good challenges 
(tests) for theoretical studies. Based on the heavy quark 
effective theory (HQET) [19-21], we can classify them 
into two doublets, (2+,3“'') and (3''',4+), the light com¬ 
ponents of which have jP = 5/2“ and jp‘ = 7/2“, re¬ 
spectively. In this paper we shall use the method of QCD 
sum rule [22, 23] to study them, which has been success¬ 
fully applied to study the ground state (5'-wave) heavy 
meson doublet (0“,1“) [24-31], the P-wave heavy me¬ 
son doublets (0+,l+) and (l+,2+) [32-36], and the D- 
wave heavy meson doublets (1“,2“) and (2“,3“) [17]. 
In this paper we shall follow the procedures used in these 
references, and study the F-wave cs heavy meson dou¬ 
blets (2+,3+) and (3+,4+). In the calculations we shall 
take into account the ©(l/mg) corrections, where mg is 
the heavy quark mass. We note that the convergence of 
this XjrnQ expansion can be problematic because F-wave 
heavy mesons (probably) have masses significantly larger 
than the heavy quark mass. However, we still hope that 
the leading terms and the ©(I/mg) corrections could 
capture sufficiently much of the most important qualita¬ 
tive physics. We shall also carefully check this conver¬ 
gence in Sec. IV. 

This paper is organized as follows. After this Introduc¬ 
tion, we construct the F-wave cs interpolating currents 
for the heavy meson doublets (2+,3+) and (3+,4+) in 
Sec. H. These currents are then used to perform QCD 
sum rule analyses in the framework of HQET both at 
the leading order and at the ©(l/mg) order. The calcu¬ 
lations are done in Sec. H and Sec. HI, and the results 
are summarized in Sec. IV. 


II. THE SUM RULES AT THE LEADING 
ORDER (IN THE mg->■ oo LIMIT) 

The heavy meson interpolating currents have been sys¬ 
tematically constructed in Refs. [32-34]. Here we fol¬ 


low Ref. [17] and briefly show how we construct the F- 
wave interpolating currents. We denote them as i 

where j and P are the total angular momentum and par¬ 
ity of the heavy meson, and ji is the total angular mo¬ 
mentum of the light components (containing three orbital 
excitations). We have the following relation 

j = jl ® SQ , (1) 

where sq = 1/2 is the spin of the heavy quark. 

To construct the F-wave interpolating currents, we 
just need to add three derivatives to the pseudoscalar 
current hy'y^q of = 0“ and the vector current hvjfiQ 
of = 1“. By doing this, the three orbital excita¬ 
tions can be explicitly written up. We act them on the 
light (strange) quark, and the obtained field has either 
= 5/2-: 


Dr X 7 ^ 3759 , 

(2) 

Dr Dr Dr x q, 

(3) 


where D/ = — {D ■ v)v^ with — igA^. Some 

other notations are: 7/ = — hy denotes the heavy 

quark field in HQET, v is the velocity of the heavy quark, 
and denotes the transverse metric 

tensor. 

We use Eq. (2) of jp’ = 5/2“ to construct the inter¬ 
polating currents coupling to_the F-wave (2+,3“*') spin 
doublet, based on hyj^q and hyj^q: 

JIPP7/2 = ^-75 X X 7^759, (4) 

Jl 77 /T = hvir X x 7^75?. (5) 

Here x and y mean that these two currents are not pure 
2+ nor 3“*', while we can project out the two pure ones: 

477/2 = \l\hy{-ifs7vt^{vt^ 

( 6 ) 

477/7 = \l\hy7{-ifs7^rvrvri])t]q, ( 7 ) 

where Sj denotes symmetrization and subtracting the 
trace terms in the sets (ai ■ ■ - aj). We note that the 
expressions of these currents have been modified to be 
consistent with Refs. [32-34]. 

Similarly, we use Eq. (3) of jf' = 7/2“ to construct the 
interpolating currents coupling to the F-wave (3“'',4+) 
spin doublet: 

477/7 = 

(8) 

477/7"^ = {9) 
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These interpolating currents are then used to perform 
QCD sum rule analyses. As discussed in Refs. [32- 
34], we do not need to investigate all of them, but just 
choose and because the calculation 

using these two currents are a bit simpler (to be techni¬ 
cally precise, we use non-symmetrized currents to calcu¬ 
late the operator product expansion (OPE) and then do 
the “symmetrization and subtracting the trace terms”). 
Moreover, we shall fix q to be the strange quark in the fol¬ 
lowing, because we are mainly studying cs heavy mesons 
in this paper. 

We follow the procedures used in Ref. [17], and assume 
\j,P,ji) to be the heavy meson state with the quantum 
numbers j, P and ji in the mg oo limit. The relevant 
interpolating field couples to it through 


(o|j: 






( 10 ) 


where denotes the decay constant, and de¬ 

notes the transverse, traceless, and symmetric polariza¬ 
tion tensor, satisfying: 


ai-aj 

Vt 



-Pi 


S'M 


aiPl 


~otj 0j I 

■9t ] 


( 11 ) 


In this expression — q^q’l/m?, and 5) denotes 

symmetrization and subtracting the trace terms in the 
sets and (/3i---/3j). Based on Eq. (10), we 

can construct the two-point correlation function 


■■■CKj ,/3i 


Akx 


= i-^ys' 


9t '"Qt 


HjT Ji (^)) 


( 12 ) 


and calculate it at the hadron level: 

f2 . 

^J,PJi(^) = - -h higher states, (13) 

2Apj^ — LJ 

where uj = 2v ■ k denotes twice the external off-shell en¬ 
ergy. = Aj,_i/2,Pji = ^ 3 i+i/ 2 ,p, 3 i is defined to be 

Apj, = lim (pij.Pj, - mg), (14) 


where is the mass of the lowest-lying state which 

{x) couples to. 

We can also calculate Eq. (12) at the quark and gluon 
level using the method of QCD sum rule in the framework 
of the heavy quark effective theory, i.e., we insert Eq. (7) 
and (9) into Eq. (12), perform the Borel transformation, 


and then obtain (see Refs. [17, 32-35] for details): 

n3, + .5/2(Wc,T) = 4_5/2e-2A+.=/AT 

3 R 3ms 7 3ms 


2^^ L 17920^2 

(g^sGG) , 

(jJ 


LJ 


89607r2 




12807r2 


w 




10 . 


1447r2 

n4. + .7/2(Wc,T) = 4,7/2e-2A+,7/2/T 

3 R 3ms 7 3ms 


(15) 


179207r2 


19(5?GG) 4 

-UJ 


LJ 


30727r2 


89607r2 
G^dio. 


00^- 


12807r2 


U! 


(16) 


These two sum rules for (2+,3''') and (3+, 4+) are simi¬ 
lar. Similarly to Ref. [17], the quark condensate {qq) and 
the mixed condensate {gsqo'Gq) both vanish, making the 
convergence of Eqs. (15) and (16) quite good. This can 
be easily verified because we need to apply as many as 
six covariant derivatives to the light quark propagator 


zSf(y,x)^(0|T[g“(y)g-''(x)]|0) (17) 

i(5“^(y-i) 5°-'^{y-xy _ 


S’^^nio 


47r2(y — a:)2 


-{qq){y-x)' 


^ab 


EcG' 


48 

1 


327r2 2 '^"(2 /-x)2 


8772 (y — a:)2 

{a^^y-x) + (y-a:)a^’^) 


1 


ab 


ScG 


47r2 2 ^‘'{y-xy 


{y - x)yfj.x^ 


Differently, we need to carefully deal with the gluon terms 
contained in these covariant derivatives in order to eval¬ 
uate the gluon condensate, which gives significant con¬ 
tribution. The gluon condensate and the strange quark 
mass take the following values [32-35, 37]: 


(y^GG) = 0.005 ± 0.004 GeV^ (18) 

TT 

m, = 0.15 GeV. (19) 


We note that the radiative corrections are not taken 
into account in our calculations, which can be impor¬ 
tant but not easy to evaluate, because the six covariant 
derivatives also contribute to them (see discussions re¬ 
lated to /p in Ref. [38] and related references). However, 
we expect that they would lead an uncertainty signifi¬ 
cantly smaller than the gluon condensate and the charm 
quark mass. Hence, we shall discuss the change of the 
latter two parameters in Sec. IV, but do not discuss the 
radiative corrections any more. 

To obtain Apj-,, we just need to differentiate 
Log[Eq. (15)] and Log[Eq. (16)] with respect to [—2/r]: 


Apj, {uJc,T) 


^[Lognj.pj, juJcT)] 
d[-2/T] 


( 20 ) 
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Then we can use it to further evaluate fpj ^: 

fpM^c.T) = V^e2Ap..,(-<=T)/Txn,,p,,,K,T).(21) 

There are two free parameters in these equations: the 
Borel mass T and the threshold value uJc- We need to 
fix these two parameters to evaluate Apj,{uJc,T) and 
fpji 



FIG. 1: The variations of CVG and PG with respect to the 
Borel mass T. The sum rule (15) for the current ^ is 

used in both figures. 


We use two criteria to fix the Borel mass T. One cri¬ 
terion is to require that the high-order power corrections 
be less than 30% to determine its lower limit 


Convergence (CVG) = 


Tjhigh—order 


(oo,r) 


nj.pji (oo, T) 


< 30%, 


( 22 ) 

where T) denotes the high-order power 

corrections, for example, 


Tjhigh—order 
^^ 3 ,+ , 5/2 


+C+) 



\ ialGG) 4 I 

J 2ms 

1447r^ 




(23) 


The other criterion is to require that the pole contribu¬ 
tion (PC) be larger than 10% to determine its upper limit 


pC=B2£:2i(+l^>10%. 


(24) 


Altogether we obtain a Borel window Tmin < T < T^ax 
for a fixed threshold value uJc- This uJc is the other free 
parameter, which will be fixed in Sec. IV. Here we pro¬ 
ceed using the sum rule (15) and taking ujc = 3.0 GeV 
as an example. Using this value of Wc, we obtain a Borel 
window 0.376 GeV < T < 0.513 GeV for the sum rule 
(15): the lower limit is determined by using the first crite¬ 
rion of CVG, as shown in the top panel of Fig. 1, and the 
upper limit is determined by using the second criterion 
of PC, as shown in the bottom panel of Fig. 1. 

Finally, we show the variations of A_|__ 5/2 f+^ 5/2 

with respect to the Borel mass T in Fig. 2. We show them 
in a broader region 0.3 GeV < T < 0.6 GeV, while these 
curves are more stable in the Borel window 0.376 GeV 
< T < 0.513 GeV. We obtain the following numerical 
results: 

A+, 5/2 = 1.40 GeV, /+, 5/2 = 0.20 GeV®/^ , ( 25 ) 

where the central value corresponds to T = 0.445 GeV 
and oJc = 3.0 GeV. 




Borel Mass [GeV] 

FIG. 2: The variations of A+, 5/2 (top) and /+, 5/2 (bottom) 
with respect to the Borel mass T. In both figures we take oJc 
= 3.0 GeV and the Borel window is 0.376 GeV < T < 0.513 
GeV. 


The procedures are the same for different values of uJc- 
We give it a large range 2.5 GeV< oJc < 3.5 GeV, but find 
that there are Borel windows as long as sq > 2.7 GeV^. 
The corresponding Borel windows and the numerical re¬ 
sults of A+, 5/2 and +, 5/2 are listed in Table. I. We note 
that this table is shown in Sec. IV, where we shall fix ojc 
to evaluate m, 2 ,+, 5/2 and TO 3 ,+, 5 / 2 . 
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Similarly, we use the sum rule (16) to perform QCD 
sum rule analyses. The Borel windows and the numerical 
results of A_|_ 7/2 and f+j /2 for various values of Wc are 
listed in Table II, also shown in Sec. IV. Here we show 
the variations of A_|_ 7/2 and f+j /2 with respect to the 
Borel mass T in Fig. 3, when we take ujc = 3.0 GeV and 
the Borel window is obtained to be 0.365 GeV < T < 
0.518 GeV. Again these curves are more stable inside 
this window. We obtain the following numerical results: 

A+, 7/2 = 1.37 GeV, /+, 7/2 = 0.19 GeV®/^ , (26) 

where the central value corresponds to T = 0.442 GeV 
and uJc = 3.0 GeV. 



Borel Mass [GeV] 

FIG. 3: The variations of A+_ 7/2 (top) and f+, 7/2 (bottom) 
with respect to the Borel mass T. In both figures we take uj^ 
= 3.0 GeV and the Borel window is 0.365 GeV < T < 0.518 
GeV. 


III. THE SUM RULES AT THE ©(l/mg) ORDER 

In the previous section we have calculated Apj, = 
limmQ-).oo(TOj_pj, — mg), the value of which is the same 
for both Ajj_ 7 / 2 .fj'i s-ad Ajjp 7 / 2 .f • To differentiate the 
masses within the same doublet, i.e., between 
and m,jjpi/ 2 ,pj-,, we need to work at the 0{l/mQ) or¬ 
der, which will be done in this section. Again we follow 
the procedures used in Ref. [17] (see Refs. [17, 32-35] 
for details), and write the pole term on the hadron side, 


Eq. (13), as: 

if + Sf? 

^ _P 25mp 2J6f 
2A-Lu (2A - a;)2 2A - w ’ ^ 

where we have omitted the subscripts j, P, ji for simplic¬ 
ity. The corrections to the mass can be evaluated 

through 

Srrij^Pj^ = —-^--—(Kpji+djj^Cmag'P‘P,ji) , (28) 

where dj,- 1 / 2 , 7 , = + 2 , d/,+ 1 / 2 , 7 , = and 

Cmag{jnQ/^j) = [og( toq)/ as( m)]^'^^° with /3o = 11 - 
2ny/3. The two corrections Kp,ji and Spj-, come from 
the nonrelativistic kinetic energy and the chromomag- 
netic interaction, respectively. We can calculate them 
using the method of QCD sum rule in the framework of 
HQET. We obtain the following two equations for df+_ 5/2 
and E_|_ 5 / 2 : 


jp -2IV^572jT 

/+,5/2^-l-,5/2e 

1 70 161(g,^GG) ^ 

92167r2 307207r2 

/+,5/2S+,5/2e-“+-^/^/^ = 


'2ms - 


pUJc 

\{9lGG) e] 

J 2ms 

- w 

25600772 


e-‘^/'^duj , (29) 

(30) 


and the following two equations for iF+ 7/2 and E+ 7 / 2 : 

/+,7/2^+,7/2e-“+'VVT 


f 


3 ^ 2\l{g^GG) 


358407r2 


/+,7/2^-l-,7/2e 


-2A 


+ . 7/2 


/T _ 


61440772 

{alGG) . 


26800772 


w 


e . 

(32) 


Again, these sum rules for (2+, 3'*') and (3+, 4+) are sim¬ 
ilar. Then A^+, 5 / 2 , ^+, 5 / 2 ? dCp 7 / 2 , and E_|_ 7^2 can be 
simply obtained by dividing these equations with respect 
to the sum rules (15) and (16). We evaluate their numer¬ 
ical results in the Borel windows derived in the previous 
section, and list them for various values of uic in Tables I 
and II. 

Here we take uic = 3.0 GeV as an example, and show 
their variations with respect to the Borel mass T in 
Figs. 4 and 5. We use the Borel windows 0.376 GeV 
<T < 0.513 GeV for ATp 5/2 and S_,__ 5 / 2 , and obtain the 
following numerical results: 

K+,5/2 = -4.23 GeV^ , S+, 5/2 = 0.014 GeV" , (33) 

where the central value corresponds to T = 0.445 GeV 
and (jJc = 3.0 GeV. We use the same Borel window 0.365 
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FIG. 4: The variation of i('+_ 5/2 (top) and E+, 5/2 (bottom) 
with respect to the Borel mass T. In both figures we take cUc 
= 3.0 GeV and the Borel window is 0.376 GeV < T < 0.513 
GeV. 


GeV <T < 0.518 GeV for i^+, 7/2 and E+^ 7 / 2 , and obtain 
the following numerical resnlts: 

K+^t/2 = -3.25 GeV^, S+^ 7/2 = 0.012 GeV^ , (34) 

where the central value corresponds to T = 0.442 GeV 
and uJc = 3.0 GeV. 


IV. NUMERICAL RESULTS AND 
DISCUSSIONS 

The mass of the F-wave cs heavy mesons can be ob¬ 
tained using Eqs. (14) and (28). We use (mu*^, mu^j) 
to denote the mass of the heavy mesons belonging to the 
(2+,3+) spin doublet, and they satisfy: 


1 ^ 


1 

—( 577777.3 + 777777,3) 

— 777c + A_|__5/2 — 

4777c 



(35) 


3 y 

— - ^+,b/2- 

777c 

(36) 


We use (tod' 3 , "m ) to denote the mass of the heavy 
mesons belonging to the (3+,4+) spin doublet, and they 




Borel Mass [GeV] 

FIG. 5: The variation of (top) and (bottom) 

with respect to the Borel mass T. In both figures we take cJc 
= 3.0 GeV and the Borel window is 0.365 GeV < T < 0.518 
GeV. 


satisfy: 

^(777777-3+977777.4) = mc + A+,7/2 - ^A+, 7 / 2 , 

(37) 

4 

77777. - m 77 ' = —7/2 . ( 38 ) 

s 4 S 3 

In this paper we use the charm quark mass m,. = 1.275 ± 
0.025 GeV, which is evaluated in the MS scheme [1]. 
Using the above equations, we calculate ( 77777 .^, 7777743 ), 
(77777-3,77777.4), their differences for various threshold 

values ijJc- The results are listed in Tables I and II. For 
completeness, we also list Borel windows, fp, 31 1 

Kp j ^, and Epj-, for various Wc. 

Now we can fix the the threshold value uJc- Our cri¬ 
terion is to require that the Wc dependence of the mass 
prediction be the weakest. We show variations of 77777 *^ 
and 77777.4 with respect to the threshold value Wc in the 
top panels of Figs. 6 and 7, and quickly notice that 
this dependence is the weakest around ojc ^ 2.8 GeV 
for both cases. Accordingly, we choose the region 2.7 
GeV< UJc < 3.0 GeV as our working region. We ob¬ 
tain the following numerical results for the (2’'', 3'*") spin 
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FIG. 6: The variations of (top) and Am_|__ 5/2 (bottom) 

with respect to the threshold values u)c- The upper and lower 
bands are obtained by using Tmin and Tmax, respectively. 



FIG. 7: The variations of (top) and Am+_ 7/2 (bottom) 

with respect to the threshold values cJc. The upper and lower 
bands are obtained by using Tmin and Tmax, respectively. 


doublet: 

= 3.45 ± 0.25 GeV, 

ruD^s = 3.50 ±0.26 GeV, (39) 

Am+, 5/2 = 0.046 ±0.030 GeV, 

where the central value corresponds to T = 0.418 GeV 
and ujc = 2.8 GeV. Here the uncertainties are due to the 
Borel mass T, the threshold value Wc, and the uncertainty 
of the gluon condensate {^GG) = 0.005 ± 0.001 GeV"^. 
We obtain the following numerical results for the (3“'', 4+) 
spin doublet 

= 3.20±0.22 GeV, 

niD-^ = 3.26±0.23 GeV, (40) 

Am+, 7/2 = 0.053 ±0.044 GeV, 

where the central value corresponds to T = 0.417 GeV 
and oJc = 2.8 GeV. However, we note that the mass differ¬ 
ences within the same doublets, Am_|_ 5/2 = 
and Am_|_ 7/2 = ~ depend on the thresh¬ 

old value Wc, as shown in the bottom panels of Figs. 6 
and 7. 



FIG. 8: The variations of with respect to the threshold 
values uJc- The solid curve is our previous result which has 
been plotted in the top panel of Fig. 6; the short-dashed curve 
is obtained by using the gluon condensate {^GG) = 0.012 
GeV* [22]; the long-dashed curve is obtained by using the 
pole mass of the charm quark me = 1.67 GeV [1]. 


The above analyses suggest that there is a heavy me¬ 
son spin doublet (2+,3+) whose masses are around 3.45 
GeV and 3.50 GeV, and a spin doublet (3+,4+) whose 
masses are around 3.20 GeV and 3.26 GeV. The latter 
is consistent with recent theoretical studies [ 6 , 12, 13], 
while the former is larger but still within uncertainties, 
as shown in Table HI. We note that the two sum rules for 
(2+,3+) and (3''',4''") are similar, see Eqs. (15) and (16), 
Eqs. (29) and (31), and Eqs. (30) and (32), so the mass 
difference between (2+, 3+) and (3+, 4+) may be (partly) 
due to the theoretical uncertainty of the numerical anal¬ 
ysis. Moreover, the expansion on the charm quark mass 
for the (2+,3“'') spin doublet is 

Eq. (35) ~ me ± 1.4 GeV ± 1.0 GeV, 


(41) 




















TABLE I: The mass of the heavy mesons belonging to the (2^,3”'') spin doublet and and their differences 

Am+, 5/2 = moss ~ for various threshold values tOc. We also list Borel windows, A+, 5 / 2 , f+,5/2, K+^5/2, and E+^ 5/2 for 

completeness. 


cjc[GeV] Borel window[GeV] A[GeV] /[GeV®/^] A[GeV®] 

E[GeV®] 

mu., [GeV] 

mu„ 3 [GeV] Am[GeV] 

2.7 

[0.376,0.426] 

1.4835 

0.1748 

-3.9886 

0.02557 

3.5055 

3.5656 

0.0601 

2.8 

[0.376,0.459] 

1.4210 

0.1701 

-3.9779 

0.01961 

3.4491 

3.4953 

0.0462 

2.9 

[0.376,0.488] 

1.4028 

0.1812 

-4.0794 

0.01613 

3.4556 

3.4935 

0.0379 

3.0 

[0.376,0.513] 

1.4036 

0.1992 

-4.2317 

0.01375 

3.4895 

3.5218 

0.0323 

3.1 

[0.376,0.537] 

1.4151 

0.2221 

-4.4156 

0.01203 

3.5394 

3.5677 

0.0283 

3.2 

[0.376,0.560] 

1.4333 

0.2491 

-4.6213 

0.01071 

3.5997 

3.6250 

0.0253 

3.3 

[0.376,0.582] 

1.4555 

0.2800 

-4.8432 

0.009658 

3.6669 

3.6897 

0.0228 

3.4 

[0.376,0.603] 

1.4807 

0.3146 

-5.0789 

0.008789 

3.7395 

3.7602 

0.0207 

3.5 

[0.376,0.624] 

1.5081 

0.3532 

-5.3259 

0.008076 

3.8163 

3.8353 

0.0190 

TABLE 11: The mass of the heavy mesons belonging to the (3^,4^) spin doublet and 

for various threshold values u!c- We also list Borel windows, A_|_ 7/2, /+,7/2i 

completeness. 

, and their < 
K+,7/2, and 

cjc[GeV] Borel window[GeV] A[GeV] /[GeV®^] A"[GeV=] 

E[GeV®] 

mu '3 [GeV] 

mu* 4 [GeV] Am[GeV] 

2.6 

[0.365,0.404] 

1.4690 

0.1594 

-3.0566 

0.02796 

3.2940 

3.3817 

0.0877 

2.7 

[0.365,0.440] 

1.3884 

0.1486 

-2.9815 

0.02076 

3.2114 

3.2765 

0.0651 

2.8 

[0.365,0.469] 

1.3649 

0.1571 

-3.0294 

0.01682 

3.2043 

3.2570 

0.0527 

2.9 

[0.365,0.495] 

1.3632 

0.1725 

-3.1262 

0.01420 

3.2262 

3.2707 

0.0445 

3.0 

[0.365,0.518] 

1.3735 

0.1926 

-3.2523 

0.01234 

3.2645 

3.3032 

0.0387 

3.1 

[0.365,0.542] 

1.3907 

0.2165 

-3.3976 

0.01090 

3.3126 

3.3469 

0.0343 

3.2 

[0.365,0.564] 

1.4127 

0.2442 

-3.5579 

0.009792 

3.3680 

3.3988 

0.0308 

3.3 

[0.365,0.585] 

1.4379 

0.2754 

-3.7294 

0.008901 

3.4284 

3.4563 

0.0279 

3.4 

[0.365,0.606] 

1.4651 

0.3102 

-3.9120 

0.008126 

3.4928 

3.5183 

0.0255 

3.5 

[0.365,0.627] 

1.4941 

0.3488 

-4.1013 

0.007510 

3.5600 

3.5836 

0.0236 

while the expansion 

for the (3+,4+) 

spin 1 

doublet has 

Ref. [37, 39] for detailed discussions). 


better convergence 

Eq. (37) - TOc + 1.4 GeV + 0.8 GeV, 


(42) 


This suggests that our results for the latter doublet are 
more reliable. 

To make our analyses complete, we try to change the 
values of the parameters used in the previous analyses 
and redo the calculations: 

1 . As shown in sum rules (15) and (16), the gluon con¬ 
densate is important. Besides the value listed in 
Eqs. (19), {^GG) = 0.005 ±0.004 GeV^ [37, 39], 
the value {^GG) = 0.012 ± 0.004 GeV^ is also 
widely used in QGD sum rule studies [22] (see 


We use 

this value and redo the numerical analyses. The 
mass of the 2~^ heavy meson, is shown in 

Fig. 8 with respective to the threshold value Wc, 
using short-dashed curves. The obtained result is 
even larger than 4.0 GeV, which is not very reli¬ 
able / reasonable. 

2. We change the charm quark mass from the MS 
value rric = 1.275 ± 0.025 GeV to its pole mass 
rric = 1.67 ± 0.07 [1], and redo the numerical anal¬ 
yses. The result is shown in Fig. 8 with respective 
to oJc using long-dashed curve. The obtained result 
is about 200 MeV larger than our previous result, 
suggesting that our results for the masses of the 
heavy mesons can have signihcant theoretical un- 
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TABLE III: Masses of _F-wave charmed-strange (cs) mesons (in GeV). 


State This work 

State Ref. [6] Ref. [12] Ref. [13] 

2+ in (2+,3+) 3.45 ±0.25 

UF 2 

3.159 

3.230 

3.224 

3+ in (2+,3+) 3.50 ±0.26 

l^Fg 

3.151 

3.266 

3.247 

3+ in (3+,4+) 3.20 ±0.22 

l^Fg 

3.157 

3.254 

3.203 

4+ in (3+,4+) 3.26 ±0.23 


3.143 

3.300 

3.220 


certainties (see also discussions in Ref. [17]). 

We can similarly replace the charm quark by bot¬ 
tom quark and study the bs system (the factor Cmag in 
Eq. (28) is taken to be 0.8 [34, 35]). Again, these masses 
depend much on the bottom quark mass mb, whose value 
has large uncertainties. When we use the IS mass value 
mb = 4.66 GeV [ 1 ], we can obtain the mass of the F-wave 
bs heavy-light mesons to be around 6.3 GeV, consistent 
with the results obtained in Ref. [12, 14]. Their mass dif¬ 
ferences are 0.010 GeV and 0.014 

GeV. However, if we replace the strange quark by up and 
down quarks, the sum rules (15) and (16) would become 
too simple to investigate the non-strange F-wave heavy 
mesons. 

In summary, in this work we adopt the QSR ap¬ 
proach to study the mass spectrum of F-wave heavy-light 
mesons in the framework of HQET. We obtain two simi¬ 
lar sum rules for (2+,3“'') and (3+,4+), see Eqs. (15) and 
(16), Eqs. (29) and (31), and Eqs. (30) and (32). Our 
results suggest that there is a cs heavy meson spin dou¬ 
blet (2+,3+) whose masses are 3 +^ = (3.45 ± 0.25, 

3.50±0.26) GeV, with mass difference Am[j]® 3^2 = 0.046± 
0.030 GeV, and a spin doublet (3+, 4+) whose masses are 
mf^i 4 +) = (3.20 ±0.22, 3.26 ±0.23) GeV, with mass dif¬ 
ference Am[^®j ^2 = 0.053 ±0.044 GeV. We note that this 
is a pioneering work and these results are provisional. 

Finally, we would like to note that this is a pioneer¬ 
ing study applying HQET-based QSR to study F-wave 
heavy-light mesons (see also discussions in Sec. I). They 


have large orbital excitations F = 3, which can be explic¬ 
itly written as three covariant derivatives, and are not 
easy to deal with. However, because the LHCb experi¬ 
ments have just observed F-wave heavy mesons [3, 4], the 
theoretical analyses on F-wave heavy mesons, including 
our study in current paper, become helpful to the fur¬ 
ther experimental exploration of them. Moreover, they 
are also good challenges (tests) for theoretical studies. 
In the following experiments such as LHCb and Bellell, 
searching for higher excitations of heavy-light mesons will 
be an important task. We expect more experimental and 
theoretical progresses on higher excitations of heavy-light 
mesons, which will make our knowledge of heavy-light 
meson family become more and more abundant. This 
will improve our understanding to the non-perturbative 
behavior of QCD, and inspire ideas for the improvement 
of QCD sum rule itself. 
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